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Abstract: Playing an important role in global warming and plant growth, relative humidity (RH) has 

profound impacts on production and living, and can be used as an integrated indicator for evaluating the 

wet-dry conditions in the arid and semi-arid area. However, information on the spatial-temporal variation 

and the influencing factors of RH in these regions is still limited. This study attempted to use daily 

meteorological data during 1966—2017 to reveal the spatial-temporal characteristics of RH in the arid region 

of Northwest China through rotated empirical orthogonal function and statistical analysis method, and the 

path analysis was used to clarify the impact of temperature (T), precipitation (P), actual evapotranspiration 

(ET), wind speed (W) and sunshine duration (S) on RH. The results demonstrated that climatic conditions 

in North Xinjiang (NXJ) was more humid than those in Hexi Corridor (HXC) and South Xinjiang (SXJ). 

RH had a less significant downtrend in NXJ than that in HXC, but an increasingly rising trend was observed 

in SXJ during the last five decades, implying that HXC and NXJ were under the process of droughts, while 

SXJ was getting wetter. There was a turning point for the trend of RH in Xinjiang, which occurred in 2000. 

Path analysis indicated that RH was negatively correlated to T, ETa, W and S, but it increased with increase 

of P. S, T and W had the greatest direct effects on RH in HXC, NXJ and SN), respectively. ET, was the 

factor which had the greatest indirect effect on RH in HXC and NXJ, while T was the dominant factor in 

SX]. 
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1 Introduction 


Water vapor is one of the most vital components in the atmosphere for it is not only a necessary 
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condition for the formation of precipitation, but also an important greenhouse gas, which is a key 
driving force of climate change (Held and Soden, 2006; Dessler et al., 2008; Li and Zha, 2018). 
There were many climatic prediction models possessing strong water vapor feedback (Held and 
Soden, 2000). However, most of them did not give an attention to water vapor component in the 
models, as most researchers attribute changes of water vapor to its natural dynamics. Previous 
studies mainly focused on dew-point temperature (Robinson, 1998, 2000), specific humidity 
(Willet et al., 2007; Iwasaki, 2016) and wet-bulb temperature (Moratiel et al., 2017; Kang and 
Eltahir, 2018). In contrast, relative humidity (RH), as a parameter from direct observation, can 
reflect the combined effect of other meteorological variables (e.g., precipitation and temperature) 
which may be used to assess regional wet-dry conditions (Gaffen and Ross, 1999; Lu and Takle, 
2010). 

Only a few studies focused on RH over China, even all over the world (Li and Zha, 2018). 
Although there were different trends of RH among different regions, it kept globally unchanged 
(Dessler et al., 2018). Dai (2006) opined that upon the application of grid dataset, the global trend 
of near-surface RH was insignificant, but it showed a remarkably positive trend in western China, 
United States and India. It is similar to the finding from Willett et al. (2008) that changes in RH 
were not significant except for some seasonal changes over the globe, tropics and Northern 
Hemisphere during 1973-2003. Surratt et al. (2004) and Vincent et al. (2007) systematically 
analyzed long-term tendency of RH over the United States and Canada. Results indicated that RH 
showed a downward trend in some areas, but overall, RH was significantly rising over the North 
America. A significant downward trend of RH was observed in North China, Northeast China and 
the eastern Northwest China. Likewise, a notable upward trend was eminent in the Qinghai-Tibet 
Plateau and the western Northwest China during the last 50 years (Song et al., 2012). 

Researches on the variations of RH in lower atmosphere are essential to eco-geographical 
environment, as RH is closely related to the hydrological circulations (Gedney et al., 2006) and the 
surface energy exchanges (Xie et al., 2011). The significant increase of temperature in arid and 
semi-arid areas, especially in semi-arid regions, contributes about 44% of the global land-surface 
warming (Huang et al., 2012). Trenberth et al. (2003) indicated that atmospheric moisture is 
expecting to rise with temperature, which is the same as the research that lower-tropospheric 
temperature has increased with RH holding steady in the United States (Wang et al., 2008). With 
climatic warming, the amount of moisture in the atmosphere is expected to increase much faster 
than the total precipitation amount (Xie et al., 2011). Therefore, studies on climate change in arid 
and semi-arid regions are the non-negligible part of the global change research. RH, as a 
comprehensively meteorological factor, plays a critical role in climate change in the arid region of 
Northwest China (ANC). However, the dynamics of near-surface RH over ANC under climatic 
change and its influencing factors remain unclear. The direct effect and indirect effect of 
precipitation, temperature, wind speed and sunshine duration on RH are still not clearly clarified. 
It focused on one-dimensional or multi-dimensional statistical methods and lacked analyses about 
the impact path in previous researches. The direct effect of a factor on RH may be positive, but it 
may turn into a negative effect through other factors. To further explore the interactive relationship 
of climatic parameters with RH, path analysis method was used in this study to investigate the role 
of the climatic factors on RH in the study area. The findings of the study are able to provide 
evidences for water resources management and agricultural practices. 

The objectives of this study were as follows: (1) to reveal the spatial-temporal characteristics of 
RH over ANC by using rotated empirical orthogonal functions (REOF) from 1966 to 2017; (2) to 
identify intra-annual variation of RH in sub-region scale through heat map; and (3) to clarify the 
impacts of meteorological variables on RH in each sub-region by path analysis. Some different 
spatial distributions and temporal characteristics between these sub-regions result from REOF 
analysis may be obtained from these analyses. Except that, it can get how temperature, precipitation 
and wind speed affect RH, positive or negative effect? There are some new research conclusions 
which will lay a scientific foundation for revealing the variations of atmospheric water vapor 
content under the global climate change. And they have certain scientific guiding significance for 
regional water resources management and spatial layout of ecological construction. 
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2 Study area and data 


2.1 Study area 


This study takes ANC as the study area (34°-50°N, 72°-108°E), including Xinjiang, the western 
Inner Mongolia, northwestern Gansu and part of Ningxia (Fig. 1). It covers an area of 2.9x10° km?, 
accounting for about 27% of total Chinese terrestrial area. As one of the driest areas in the world 
(Yao et al., 2018), ANC has a typical continental climate with low annual precipitation (less than 
200.0 mm) (Chen et al., 2014; Deng et al., 2014) and a wide range of annual mean temperature in 
different area (from —2°C to 19°C) (Liu et al., 2010). The precipitation exhibited an increasing trend 
in ANC in recent years, with a rate of 7.7 mm/10a. The maximum precipitation occurred in 2016 
(200.2 mm), while the minimum in 1997 (99.8 mm). Complex topographic components including 
desert, forest, farmland, grassland and mountain glaciers. There are three main sources of water 
vapor, which are Arctic Ocean, Atlantic Ocean and Indian Ocean (Wang and Ren, 2005). The 'V'- 
shaped terrain towards the west catches more moisture content in air from Atlantic Ocean through 
westerly jet, which forms rainfall or snowfall on the windward slope (Guo and Li, 2015). ANC is 
in the transitional zone between humidification in eastern Eurasian continent and drying in western 
continent, which is sensitive to the climate change (IPCC, 2007). The Tianshan Mountains, Altay 
Mountains, Kunlun Mountains and the north Qilian Mountains are located in the study area. The 
Tarim Basin, Junggar Basin and Hexi Corridor (HXC) are distributed among these mountains. 
Tianshan Mountains plays an important role in natural conditions in this area. Precipitation in the 
North Xinjiang (NXJ; north of Tianshan Mountains) is greater than that in the South Xinjiang (SXJ; 
south of Tianshan Mountains) (Guo and Li, 2015). 
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Fig.1 Distribution of meteorological stations in the arid region of Northwest China 


2.2 Data 


Daily meteorological records were selected for a 52-a period from January 1966 to December 2017 
in this study. Collected data comprised daily mean RH, mean temperature , mean wind speed at 10 
m height, precipitation and sunshine duration, as well as maximum temperature, minimum 
temperature, mean atmospheric pressure, which were provided by National Meteorological 
Information Center (http://data.cma.cn/). The data consistency was checked by the double mass 
curve method (Zhang et al., 2012). There are more than 100 meteorological stations in ANC, but 
some of them exist missing data. Stations with the missing data sets which were over 5% or more 
than 5 d in a month were removed in this study. After this process, only 81 stations (Tables S1-53) 
are available in this study (Fig. 1). In other cases, missing data were interpolated by their 
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neighboring days or the neighboring stations of the same dates. Daily actual evapotranspiration 
(ETa) was calculated by Advection-Aridity (AA) model. The monthly and annual data were then 
obtained by these variables. 


3 Methods 
3.1 AA model 


AA model, which was based on the complementary theory, was applied to estimate the actual 
evapotranspiration in this study. In the AA model, potential evapotranspiration was calculated 
through the Penman method (Penman, 1948), and the evapotranspiration in wet environment was 
obtained by the Priestley-Taylor equation (Priestley and Taylor, 1972; Brutsaert and Stricker, 1979). 
The equation (Su et al., 2017; Jian et al., 2018) for AA model was as follows: 


k k 
ET, =(2a Do G) mo (1) 
where ET, is the daily actual evapotranspiration (mm); a is the Priestley-Taylor evaporation 
coefficient that varies spatially in relation to the characteristic of the underlying surface; k is the 
slope of the saturation vapor pressure curve at the air temperature (kPa/°C); y is the psychometric 
constant; Rn is the net radiation near the surface expressed as the equivalent vaporization rate 
(MJ/(m?-d)); G is the soil heat flux (mm/d); and Ea is the drying power of the air (mm/d) which 
represents the capacity of the ambient air to extract water from the surface and is an increasing 
function of the vapor pressure deficit of the air and of wind speed. 
AA model can achieve highly accurate results only based on relatively few parameters and a 
simple algorithm. It has been widely used in different climatic regions, including arid and semi- 
arid regions (Haque, 2003; Liu et al., 2010). 


3.2 Rotated empirical orthogonal function (REOF) 


In statistics, the empirical orthogonal function (EOF) analysis is known as the principal component 
analysis, but the EOF method can find both time and spatial patterns. The first a few eigenvectors 
of EOF can represent the variability structure of the fields of climatic variable (Lorenz, 1956). 


Xmxn=EOFmx px PCpxn, (2) 


where X is the matrix of empirical orthogonal function; m is the number of points in the study area; 
nis the number of days which starts from 1966 to 2017; EOF is the matrix of the spatial pattern; 
PC is the matrix of temporal coefficients; and p is the number of spatial patterns or temporal 
coefficients of EOF analysis. 

However, it is difficult to interpret the spatial characteristics of the different regions through EOF 
analysis. Thus, REOF was introduced to circumvent these problems (Richman, 1986). The most 
popular type of REOF is the varimax REOF analysis (Kaiser, 1958). Spatial field decomposed by 
REOF would concentrate the high values of loading factors in a small region, while loading factors 
in most of regions are close to zero. As a result, it will make spatial structure clear and can be used 
for geographical zoning. Based on results from its strong spatial interpretation, REOF has been 
extensively used in atmospheric, oceanic and climatic researches (Lian and Chen, 2012; Wang et 
al., 2014, 2015). More information on REOF analysis and its detailed algorithm can be found in 
Hannachi et al. (2007). 

In this study, after rotating the first 11 spatial modes and their time functions which variance 
contributions occupy more than 85% of the total variance contributions, it can get the rotated 
component loading patterns (REOFs) and rotated principal components (RPCs), which reflect 
spatial pattern and time evolution of a variable, respectively. Then the first four REOFs, which 
occupied about 64.81% of the total variance, were shown by the inverse distance weighting 
interpolation with GIS technique, and sub-regions were then divided based on the spatial 
characteristics. Only REOF spatial weights above 0.30—0.35 could be interpreted as remarkable 
signals and then used to reveal the spatial features (Richman and Gong, 1999). 
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3.3 Path analysis 


Path analysis, as an extension of multi-regression method, was introduced to differentiate between 
correlation and causation, and clarify the direct and indirect effects of independent variables on 
dependent variable (Wright, 1934). Initially, path analysis focused on the field of disease of 
multifactorial inheritance, but it was extensively applied in geographical and ecological studies for 
the past few years (Li et al., 2011; Zhang et al., 2014). The structure of data could be comprehended 
through path analysis (Kozak and Kang, 2006), thus it can be used to assess effects of 
meteorological factors on RH in this study. 

Path coefficients were defined as a version of standardized linear regression weights by Wright 
(1921), which can be applied to verify causation between causal factors of structural equation 
model. Assumed that there are n independent variables, x1, X2, ..., Xn, and a dependent variable y, 
correlation coefficient between x; and y is equal to the sum of the direct coefficient and indirect 
coefficient of x; to y: 


n 
yi = Pyi a aa, 4 (3) 


where ry; and rj are the correlation coefficients between x; and y, and between x; and aj, respectively; 
and pyi is the direct coefficient of x; to y. Meanwhile, coefficient of residual term was calculated to 
indicate if the effect of these factors is strong enough. 


De = J pa Pyy > (5) 


where p is the residual path coefficient; D? is the coefficients of determination; p is the path 
coefficient; and r is the correlation coefficient. Additional details about this method are described 
in Wright (1934). The method used here for assessing the influencing meteorological factors of 
RH, which was based on daily meteorological records, can help in getting stronger appreciation of 
the internal relationships of selected variables, especially their indirect connections. 


4 Results 


4.1 Spatial and temporal characteristics of RH using REOF 


It required a dimensionality reduction through EOF analysis before applying REOF. The first 11 
EOF modes, which occupied more than 85% of the total variance contribution, were selected as 
inputs for the REOF (Table 1). After the orthogonal varimax rotations of the loading factors, the 
first four modes of REOFs (REOF-1, REOF-2, REOF-3 and REOF-4; Figs. 2a—d) and RPCs (RPC- 
1, RPC-2, RPC-3 and RPC-4; Figs. 3a—d) were used to analyze the spatial and temporal patterns. 
The results showed that there were positive and negative loading factors caused by the 
decomposition of the RH field, implying that the distributions of RH were opposite phase. The 
variation trend of RH was opposite in HXC and Xinjiang (Figs. 2 and 3). 

Figure 2a showed that the greatest value area of the REOF-1 was located in HXC. Combined 
with the downward trend in RPC-1 (Fig. 3a), there was a tendency of dryness during 1966-2017 in 
HXC. However, the REOF-2 (Fig. 2b) had positive values in most parts of the study area, and the 
highest value was in SXJ, indicating a gradual increase in wetness occurred with a corresponding 
rise of RCP-2 in this sub-region (Fig. 3b). Figure 2c exhibited larger values in NXJ than those in 
other two sub-regions. Overall, there was no significant upward trend in the last five decades (Fig. 
2c). Figures 2d and 3d showed that there were high negative value and fluctuating uptrend in 
southwest Xinjiang. It was almost the same as the REOF-2 and RPC-2, which implied that these 
two regions have similar spatial and temporal patterns for RH. 

The year 2000 was detected as the turning point for RH in Xinjiang. For the period before 2000, 
the temporal coefficient of RH was obviously decreasing, but after 2000, it started increasing in 
SXJ (Fig. 3b). The temporal coefficient of RH in NXJ showed a reverse trend (Fig. 3c). It implied 
that it exhibited a wetting trend before 2000, but it got drier after 2000 in Xinjiang. 
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Table 1 Variance contributions (%) of the empirical orthogonal function (EOF) and rotated empirical orthogoanl 
function (REOF) analyses of annual average relative humidity 


EOF REOF 
Mode Variance Accumulated variance Variance Accumulated variance 
contribution (%) contribution (%) contribution (%) contribution (%) 
1 30.22 30.22 21.11 21.11 
2 15.56 45.78 15.50 36.61 
3 12.76 58.54 14.21 50.82 
4 925 67.79 13.99 64.81 
5 4.44 72.23 5.19 70.00 
6 3.64 75.87 4.42 74.42 
7 23 78.60 3.19 77.61 
8 225 80.85 mesi 80.12 
9 RID 82.97 209 82.41 
10 1.91 84.88 225 84.66 
11 1.63 86.51 1.85 86.51 
72°E 84°E 96°E 108°E 72E 84°E 96°E 108°E 
54°N ] (b) REOF-2 of RH SEN 
50°N 50°N 
46°N 46°N 
42°N 42°N 
38°N p 38°N 
Loading factors 4 Loading factors 
= High: 0.855 | 349N = High: 0.900 |34°N 
=. Low: 0.147 30°N = Low: -0.295 | 39ony 
72°E 84°E 96°E 108°E 72°E 84°E 96°E 108°E 


54°N ] (4) REOF-4 of RH 
50°N 
46°N 
42°N 


38°N 


Loading factors Loading factors 


= High: 0.172 f 34°N — High: 0.224 
= . — . 
Low:—0.888 30°N Low: —0.868 
72°E 84°E 96°E 108°E 72°E 84°E 96°E 108°E 


Fig. 2 Spatial patterns of the first four loading factors of the rotated empirical orthogonal function (REOF) of 
relative humidity (RH) 


4.2 Mean monthly variations of RH 


The study area was divided into three sub-regions resulting from the first four spatial-temporal 
patterns of REOF analysis: NXJ, including 26 stations; SXJ, including 36 stations; and HXC, 
including 19 stations. The annual average RH in each sub-region was significantly different. The 
highest value of RH took place in Ili River Valley and Altay Mountains due to the 'V'-shaped terrain 
and windward slope, respectively. RH decreased from north to south and from west to east over the 
whole ANC. 

Figure 4 exhibited mean monthly variations of RH in each sub-region. Monthly average RH was 


less than 40% between March and June (Fig. 4a) in HXC. The highest value of RH occurred in 
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Fig. 3 Temporal patterns of the first four rotated principal components (RPC) of RH. The time series is 
represented by black solid line, 5-a moving average by red solid line and the linear trend by yellow solid line. 


winter (from December to January of the next year). The NXJ had a higher RH than the other two 
sub-regions (Fig. 4b). The RH from November to March of the next year was larger than that from 
April to October and was higher than 70% in winter in NXJ. SXJ was the driest in the study area, 
especially during April and May (Fig. 4c). There were eight months in which RH was under 50% 
in a year. The lowest value of RH was centered in April and May in the study area (Fig. 4d). RH 
exhibited strong annual fluctuations, consistent with the former result which indicated that 
variations in RH has no significant trend in NXJ. The reason for the results is that NXJ was located 
on the windward slope of the Tianshan Mountains, while SXJ lied on the leeward slope and HXC 
was far from the source of water vapor. As a result, the moisture in NXJ was higher than those in 
SXJ and HXC. 


RH (%) RH (%) 
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10 70 10 70 
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Fig.4 Mean monthly variations of RH through heat-map in (a), Hexi Corridor (HXC); (b), North Xinjiang (NXJ); 
(c), South Xinjiang (SXJ); and (d), the arid region of Northwest China (ANC) during 1966-2017 
4.3 Impacts of meteorological factors on RH 


Path coefficients were calculated to analyze the impact of temperature (T), precipitation (P), actual 
evapotranspiration (ETa), wind speed (W) and sunshine duration (S) on RH. Before applying the 
path analysis to assess the impact of independent variables on dependent variable, F-test was used 
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to verify the linear significance of these variables. As shown in Tables 2-4, all the variables had 
significant correlations with RH (significant at 0.01 level). RH was significantly related with S in 
HXC and NXJ, in which correlation coefficients reached —0.53 and —0.77, respectively. For SXJ, 
RH was highly correlated with ETa. 

S and P influenced RH with direct effects of —0.39 and 0.35 in HXC (Table 2), respectively. As 
S decreased or P increased, RH increased. It can be seen that indirect path effect from all factors to 
RH existed. For direct effects, the greatest influencing factor was S. The indirect effects of ETa and 
T through S accounted for 69.4% and 68.6% of the total influence of ET, and T on RH, respectively. 
As the only factor with positive effect, direct path coefficient of P was larger than its indirect path 
coefficient. 


Table 2 Path analysis among ETa, P, S, T, W and RH in Hexi Corridor from 1966 to 2017 


Meteorological Correlation Direct path Indirect path coefficient 
factor coefficient coefficient ET, P S T w 
ET, -0.38"" 0.054" - 0.010 0.26 —0.06 —0.010 
12 0.44™ 0.35" —0.001 - 0.14 —0.02 —0.030 
S 0.53” 0.394 0.030 0.130 - 0.03 0.050 
T —0.20™ 0.084 —0.040 0.070 0.14 - 0.020 
W 0.157 0.264" —0.002 0.040 0.08 0.01 - 


Note: ™ means the correlation coefficient is significant at the 0.01 level; ** means the path coefficient is significant at the 0.01 level; - 
means no value. ETa, actual evapotranspiration; P, precipitation; S, sunshine duration; T, temperature; W, wind speed; and RH, relative 
humidity. 

For NXJ (Table 3), the direct path coefficient of T on RH was the largest (0.61), while W was 
the lowest (0.07). All factors can affect humidity through T except P. The indirect path coefficient 
from ET, through T to RH was the biggest (0.47), implying that ET, can indirectly influence RH 
by means of T. The smallest indirect path coefficient of P from ETa to RH was the lowest, which 
illustrated that precipitation had almost no impact on humidity in NXJ. This meant that 
evapotranspiration had the biggest impact on RH through S and T, which occupied 53.1% and 
62.1%, respectively. The rate of the indirect path coefficient of W though T on RH to the total 
indirect path coefficient of W was 84.2%. P was also the only factor which had positive effect on 
RH, but it was not significant. 


Table 3 Path analysis among ETa, P, S, T, W and RH in North Xinjiang from 1966 to 2017 


Meteorological Correlation Direct path Indirect path coefficient 
factor coefficient coefficient ET, P s T Ww 
ETa 0.76" 0.137 - —0.001 0.41 0.47 —0.010 
P 0.21” 0.19 0.00 - 0.14 —0.11 —0.020 
S D 0.464 0.11 —0.060 - —0.36 —0.003 
T 2076 —0.60*" 0.10 0.030 —0.28 - —0.012 
WwW 0.15" NO 0.02 0.040 —0.02 0.13 - 


Note: “™ means the correlation coefficient is significant at the 0.01 level; ** means the path coefficient is significant at the 0.01 level; 
means no value. 


The absolute values of correlation coefficients between selected variables and RH were more 
than 0.50, except P in SXJ (Table 4). W had a significant direct effect on RH, which with direct 
path coefficient of —0.45. T had the greatest indirect influence on RH with total indirect path 
coefficient of —0.46. The indirect influence of S through ETa on RH with an indirect path coefficient 
of —0.25 accounted for 48.4% of its total coefficient, while T's effect through ETa was —0.28, 
accounting for 41.6%. The total indirect path coefficient of P was negative, resulting in a weak 
correlation to RH in this sub-region. 

The residual path coefficient was significant (0.75) in HXC, which was more than that in the 
whole Xinjiang (0.57). The result indicated that the interpretation rate of selected factors in HXC 
was less than that in Xinjiang (24.7% and 43.2%, respectively), implying that other underlying 
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factors also influenced RH in this study area. Further studies are required to better understand the 
impact of meteorological factors on RH. 


Table 4 Path analysis among ETa, P, S, T, W and RH in South Xinjiang from 1966 to 2017 


Meteorological Correlation Direct path Indirect path coefficient 
factor coefficient coefficient ET, p s T Ww 
ETa 0.67" 0.314" - 0.080 0.07 0.18 0.190 
P 0.07™ 0.34" 0.07 - 0.01 —0.06 —0.150 
S —0.51™ 0.094 0.24 —0.030 - 0.13 —0.020 
in —0.66™ 0.20%" —0.28 0.100 —0.06 - —0.220 
W 0.58" -0.45*# 0.13 0.110 0.00 0.10 - 


Note: ™ means the correlation coefficient is significant at the 0.01 level; # means the path coefficient is significant at the 0.01 level; - 
means no value. 


5 Discussion 


Compared with temperature and precipitation, RH was more sensitive to local climate change 
(Gaffen and Ross, 1999; Lu and Takle, 2010). It can accurately and objectively describe the features 
of climate change in a region because it had a notable relationship with some extreme weather 
conditions like heavy rainfalls and droughts (Song et al., 2012). RH had obvious regional 
characteristics because of differences in natural conditions, especially climatic conditions, in the 
different areas. Three major sub-regions in this study, Hexi Corridor (HXC), North Xinjiang (NXJ) 
and South Xinjiang (SXJ), experienced regular fluctuations in RH. Several studies applied the 
Palmer drought severity index (PDSI) (Wang et al., 2015), standardized precipitation 
evapotranspiration index (SPEI) (Yu et al., 2014), aridity index (Huo et al., 2013) and surface 
humidity index (Ma and Fu, 2003) to evaluate the wetness/dryness conditions and their results 
showed that the wetter trend over the Northwest China was observed in the past half century. 
However, there were different variations not only in RH but also in SPEI of 3-month time scale in 
different sub-regions based on REOF analysis (Wang et al., 2017). Xinjiang exhibited a significant 
wetting trend, while HXC was increasingly dry. Another result opined that NXJ was more humid 
than HXC and SXJ, because NXJ is closer to the Arctic Ocean and Atlantic Ocean. Water vapor 
from ocean are blocked by Altay Mountains and Tianshan Mountains. Ili valley and the south slope 
of Altay Mountains are the most humid areas because they are located in the windward slope of the 
Tianshan Mountains and Altay Mountains, while SXJ and HXC lie in the leeward slope that water 
vapor cannot reach easily. The drying conditions in HXC was largely because of the influence of 
East Asian monsoon with increasingly weakening trend since the early 1980s (Xu et al., 2010). 
That was the reason why it exhibited opposite trend between Xinjiang and HXC. It was more humid 
before 2000, resulted from the enhanced westerly which brought more water vapor to Xinjiang 
(You et al., 2011). 

There were other circulation effects that influenced the dryness/wetness conditions in this study 
area, other than the prevailing westerly and East Asian monsoon. Wang et al. (2015) was of the 
view that the area of Northern Hemisphere polar vortex, in addition to the Arctic Oscillation and 
the North Atlantic Oscillation exerted certain influence on moisture in ANC. More recently, it was 
also discovered that the Pacific Decadal Oscillation exert a somewhat similar influence on drought 
evolution in Northwest China (Wang et al., 2017). Variations in RH in SXJ and HXC were also 
affected by the Tibetan Plateau High (Wang et al., 2015). Furthermore, in the drought years of 1970, 
1985 and 1997 in Xinjiang, updraft was prevalent in south of 37°N and downdraft in north of 37°N, 
while it experienced a circulation reversal in wet years of 1969, 1987 and 2002 (Xing et al., 2010). 
Liu et al. (2016) considered that there were negative relationships between PDSI and the Nifio 3.4 
index in the eastern Northwest China. It was reported that the strong El Nifio events happened in 
1987, 1998, 2003 and 2016, and it also showed higher RH during these years, indicating that El 
Nifio had an obviously increasing impact on RH in ANC. Variations in the westerly circulation, 
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South Asian High and Western Pacific Subtropical High also played a major role in other 
meteorological factors like temperature (Yuan et al., 2009; Yuan et al., 2010) and precipitation 
(Zhao et al., 2014; Li et al., 2016; Wei et al., 2017), which are of importance on the change of RH. 
Under the comprehensive influence of various circulations, the whole arid and semi-arid region of 
Northwest China is likely to get wet in the coming years (Peng and Zhou, 2017). 

A significant upward trend for precipitation and temperature was exhibited in Northwest China 
(Shi et al., 2007; Li et al., 2013; Han et al., 2016). Several studies attributed the wetter trend to the 
increasing precipitation in the Northwest China during the past 50 years (Piao et al., 2010; Deng et 
al., 2014; Li et al., 2016). However, according to a more recent study (Chen et al., 2018), 
precipitation had less contribution than evapotranspiration, temperature and wind speed to surface 
RH in the Tianshan Mountains. The increasing precipitation had little positive effects on RH, but 
temperature (Dai, 2006), wind speed and evapotranspiration had negative effects on RH from the 
path analysis. So that there was trend of humidification from increasing precipitation only, but from 
RH, it can get that it was not more humid in the large area of Northwest China. There were many 
theoretical and modelling studies that predicted that RH will remain approximately stable globally 
as the climate warms (Held and Soden, 2000; Sherwood and Meyer, 2006; Simmons et al., 2010; 
Willett et al., 2010). The drought trend occurred in HXC, which was due to the decreasing 
precipitation and increasing temperature (Yu et al., 2014). Impacts of reference evapotranspiration 
on wet/dry conditions and drought should be given more attention (Sun et al., 2016). It was found 
that notable reduction of evapotranspiration offsets the drying tendency caused by the decrease of 
precipitation in ANC, so that it exhibited a trend of humidification (Su et al., 2018). Besides these 
factors, sunshine duration and wind speed played important roles in HXC and SXJ, respectively, 
on the basis of results of path analysis in this study. 


6 Conclusions 


Based on the daily data from 81 stations in ANC over the period 1966-2017, we analyzed the 
spatial and temporal patterns of trends in RH and its influencing factors. The main findings of this 
study were summarized below. 

HXC and NXJ were observed to exhibit a drying tendency, but SXJ was getting increasingly but 
not significantly wetter. It may be caused by a limited moisture supply from the oceans where 
evaporation has been limited by sea surface temperature that have not risen with the temperature 
over land (Simmons et al., 2010). NXJ was more humid than SXJ and HXC. The reason is that 
water vapor from Atlantic Ocean and Arctic Ocean was blocked by Tianshan Mountains so that 
SXJ was drier than NXJ. For HXC, it is farther away from the source of water vapor than NXJ so 
that the water vapor can't reach HXC. The most humid season was winter, while the driest season 
is spring in the region. It exhibited more distinct wet and dry season in Xinjiang than that in HXC. 
Only precipitation could promote the growth of RH over the whole region. That is why some 
previous researchers considered that ANC exhibited tendency of humidification based on 
precipitation (Shi et al., 2007). However, RH was influenced by other meteorological factors, 
especially actual evapotranspiration, which was the most important indirect influencing factor on 
RH in HXC and NXJ, whereas temperature had the greatest indirect effect on RH in SXJ. The 
impacts of precipitation on RH was offset largely by the variation in other factors, especially wind 
speed based on the results from path analysis. 

RH is a product generated by not only climatic factors, but also land surface properties and 
human activities. Differences in properties of the underlying surface, human activities and 
increasing glacier melting caused by warming will exert influences on RH. Thus, the study on 
influencing factors of RH in this work is still uncertain, and more profound researches would be 
needful in future works. 
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Appendik 


Station number 
52267 
52323 
52418 
52424 
52436 
52446 
52447 
52495 
52533 
52546 
52576 
52652 
52661 
52674 
52679 
52681 
52787 
53502 
53602 


Station number 


Table S1 Details of meteorological stations in Hexi Corridor 
Longitude (°E) 


Station name 
Ejin Banner 
Mazongshan 
Dunhuang 
Anxi 
Yumenzhen 
Dingxin 
Jinta 
Bayan Nuru 
Jiuquan 


Gaotai 


Alagxa Right Banner 


Zhangye 
Shandan 
Yongchang 
Wuwei 
Mingin 
Wushaoling 


Jilantai 


Alagxa Left Banner 


Station name 


Latitude (°N) 


41.95 
41.80 
40.15 
40.53 
40.27 
40.30 
40.00 
40.17 
39.77 
39.37 
39.22 
38.93 
38.80 
38.23 
37.92 
38.63 
37.20 
39.78 
38.83 


Latitude (°N) 


101.07 
97.03 
94.68 
95.77 
97.03 
99.52 
98.90 

104.80 
98.48 
99.83 

101.68 

100.43 

101.08 

101.97 

102.67 

103.08 

102.87 

105.75 

105.67 


Table S2 Details of meteorological stations in North Xinjiang 


Longitude (°E) 


Elevation (m) 

940.5 
1770.4 
1139.0 
1170.9 
1526.0 
1177.4 
1270.5 
1323.9 
1477.2 
1332.2 
1510.1 
1482.7 
1764.6 
1976.9 
1531.5 
1367.5 
3045.1 
1031.8 
1561.4 


Elevation (m) 


51053 
51059 
51060 
51068 
51076 
51087 
51133 
51137 
51145 
51156 
51186 
51232 
51238 
51241 
51243 
51288 
51330 
51334 
51346 
51352 
51356 
51357 
51365 


Habahe 
Jeminay 
Burgin 
Fuhai 
Altay 
Fuyun 
Tacheng 
‘Yumin 
Emin 
Hefeng 
Qinghe 
Alashankou 
Bole 
Tuoli 
Karamay 
Beitashan 
Wenquan 
Jinghe 
Wusu 
Paotai 
Shihezi 
Shawan 


Caijiahu 


48.05 
47.45 
47.70 
47.12 
47.73 
46.98 
46.73 
46.20 
46.55 
46.82 
46.07 
45.18 
44.90 
45.93 
45.62 
45.37 
44.97 
44.57 
44.43 
44.85 
44.32 
44.33 
44.20 


86.40 
85.88 
86.87 
87.48 
88.08 
89.52 
83.00 
82.93 
83.65 
85.75 
90.38 
82.57 
82.07 
83.60 
84.85 
90.53 
81.02 
82.82 
84.07 
85.25 
86.05 
85.62 
87.53 


532.6 
968.2 
473.9 
497.0 
735.3 
860.4 
534.9 
715.7 
523.6 
1322.1 
1218.2 
369.8 
532.2 
1094.2 
450.3 
1653.7 
1357.8 
329.2 
478.7 
337.1 
444.2 
347.9 
440.5 
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Station number 


Station name 


Latitude (°N) 


Longitude (°E) 


Continued 


Elevation (m) 


51378 
51379 
51430 
51431 
51433 
51437 
51463 
51465 
51468 
51470 
51477 
51482 


Station number 
51467 
51495 
51526 
51542 
51567 
51573 
51581 
51628 
51633 
51639 
51642 
51644 
51655 
51656 
51701 
51704 
51705 
51709 
51711 
51716 
51720 
51730 
51765 
51777 
51804 
51811 
51818 
51828 
51839 
51855 
51931 
52101 


Jimsar 
Qitai 
Qapgqal 
Yining 
Nilka 
Zhaosu 
Urumqi 
Xiaoquzi 
Daxigou 
Tianchi 
Dabancheng 
Mori 


Table S3 Details of meteorological stations in South Xinjiang 


Longitude (°E) 


Station name 
Balguntay 
Shisanjianfang 
Kumux 
Bayanbulak 
Yangqi 
Turpan 
Shanshan 
Aksu 
Baicheng 
Xayar 
Luntai 
Kuqa 
Yuli 
Korla 
Tuergate 
Artux 
Wugia 
Kashi 
Akgi 
Bachu 
Kalpin 
Aral 
Tikanlik 
Ruogiang 
Taxkorgan 
Shache 
Pishan 
Hotan 
Minfeng 
Qiemo 
Yutian 
Barkol 


44.02 
44.02 
43.83 
43.95 
43.80 
43.15 
43.78 
43.48 
43.10 
43.88 
43.35 
43.83 


Latitude (°N) 


42.77 
43.22 
42.23 
43.03 
42.08 
42.95 
42.85 
41.12 
41.78 
41.23 
41.82 
41.72 
41.35 
41.73 
40.52 
39:72 
39.72 
39.48 
40.93 
39.80 
40.05 
40.55 
40.63 
39.03 
WoT 
38.43 
37.62 
37.13 
37.07 
38.15 
36.85 
43.60 


89.17 
89.57 
81.15 
81.03 
82.52 
81.13 
87.65 
87.10 
86.83 
88.12 
88.32 
90.28 


86.03 
1673) 
88.22 
84.15 
86.57 
89.23 
90.23 
80.38 
81.90 
82.78 
84.27 
82.97 
86.27 
85.82 
75.40 
76.02 
75.25 
75.715 
78.45 
78.57 
79.05 
81.27 
87.70 
88.02 
75.23 
VUAN 
78.28 
79.93 
82.72 
85.55 
81.65 
93.05 


735.9 
7995 
605.4 
662.5 
1105.3 
1851.0 
935.0 
1874.7 
3539.0 
1942.5 
1103.5 
1272.3 


Elevation (m) 


1732.4 
721.4 
922.4 

2458 
1055.3 
39.3 
398.6 

1107.1 

12292 
980.4 
982.0 

1081.9 
884.9 
899.8 

3504.4 

1298.7 

2175.7 

1385.6 

1985.1 

1116.5 

1161.8 

1012.2 
846.0 
887.7 

3090.1 

1231.2 

1375.4 

1375.0 

1409.5 

1247.2 

1422.0 

1679.4 
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Station number 


Station name 


Latitude (°N) 


Longitude (°E) 


Continued 


Elevation (m) 


52112 
52118 
52203 
52313 
51571 
51636 
51707 
51717 
51802 
51810 
51826 


Naomaohu 
Yiwu 
Hami 

Hongliuhe 

Toksun 
Xinhe 
Jiashi 

Yopurgha 
Yengisar 
Markit 
Qira 


43.75 
43.27 
42.82 
41.53 
42.80 
41.53 
39.50 
39.25 
38.93 
38.92 
37.02 


94.98 
94.70 
93.52 
94.07 
88.63 
82.62 
76.73 
76.78 
76.17 
77.63 
80.80 


479.0 
1728.6 
737.2 
1573.8 
2.0 
1014.5 
1208.6 
1206.3 
1297-5 
1177.6 
1336.5 


